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Infrared study of the charge-ordered multiferroic LuFe 2 04 
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The reflectivity of a large LuFe2 04 single crystal has been measured with the radiation field either 
perpendicular or parallel to the c axis of its rhombohedral structure, from 10 to 500K, and from 
7 to 16000 cm -1 . The transition between the two-dimensional and the three-dimensional charge 
order at Tco = 320 K is found to change dramatically the phonon spectrum in both polarizations. 
The number of the observed modes above and below Tco, according to a factor-group analysis, is 
in good agreement with a transition from the rhombohedral space group R3m to the monoclinic 
C2/m. In the sub-THz region a peak becomes evident at low temperature, whose origin is discussed 
in relation with previous experiments. 

PACS numbers: 71.30.+h, 74.25.Gz, 78.30.-j 



in 

•4— > 
CO 



I 

C 

o 
o 



> 

00 

o 

(N 

in 

o 
o 



% 



INTRODUCTION 

The multifcrroics are materials which exhibit, in a 
common range of temperatures, both ferroelectricity and 
magnetic order (ferromagnetism, ferrimagnetism or anti- 
ferromagnetism) . Much effort has been devoted to un- 
derstand the relationship between their magnetic and 
electronic degrees of freedom, [1-|3| in view of the poten- 
tial applications. Indeed, the conventional mechanisms of 
ferroelectricity and magnetism are not compatible with 
each other, and this makes the multiferroic materials 
relatively rare. One should then invoke novel mecha- 
nisms of ferroelectricity, as those associated with specific 
charge configurations in charge-ordered materials. [4| 
One of such compounds is the mixed- valence LuFe2 04 
(LFO). Its structure belongs to the rhombohedral space 
group i?3m and includes triangular double layers of Fe-0 
stacked along the c axis. A two-dimensional Fe 2+ /Fe 3+ 
charge order (CO), which forms between about 500 K and 
Tco = 320 K in those double layers .builds up in each 
Fc-0 plane a net dipole moment. |5jiZ| Below Tco, the 
CO becomes more robust and three-dimensional (3-D). 

Dielectric measurements in low electric fields show that 
the 3-D phase of LFO is ferroelectric, |jj while X-ray 
scattering indicates that the dipole moments of adjacent 
bilayers are antiparallel, [8[ so that the net polarization is 
zero. These contradictory findings can be reconciled by 
considering that the antifcrroclcctric ground state and 
the ferroelectric state may be so close in energy [8| that 
even a small electric field stabilizes the latter one. Con- 
cerning the magnetic behavior of LFO, a two-dimensional 
ferrimagnetic order is established below Tjv — 240 K |9| . 
However, a surprising decrease in the magnetic correla- 
tion length is observed below Tt, = 175 K. [lfj Finally, 
it has been found that the dielectric constant decreases 



fields, which also cause a drop in the resistivity by several 
orders of magnitude. [13j This colossal electroresistance, 



which is also observed at room temperature, makes LFO 
attractive for its potential applications. 
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sharply both under weak magnetic |ll| and electric [12| 



FIG. 1: Magnetic susceptibility of the single crystal of 
LuFe204, as measured in zero-field cooling. 



The optical response of LFO with the radiation field 
Elc has been investigated by Xu et ai., [14J together 
with other properties, between 30 meV and 6.5 eV and 
from 4 to 540 K. At zero field the authors observe two 
main peaks in the near-infrared/visible range, which is 
dominated by electronic excitations. The band at 1.1 
eV is assigned to Fe 2+ /Fe 3+ charge transfer, and that 
at about 3.5 eV to a superposition of O—p —>Fe-d and 
O— p — kLu— s charge transfer. The edge of the 1.1 eV 
peak creates the insulating optical gap at ~ 0.5 eV. A 



third feature in the mid infrared, peaked at about 0.3 eV, 
is attributed by the authors to on-site Fe 3+ excitations. 
Li et al. [15( have instead explored the sub-THz response 
of LFO by time-domain spectroscopy, observing at low T 
two collective excitations that they assigned to a central 
mode and to a soft mode of the ferroelectric phase. 

In the present experiment we have investigated both 
the sub-THz and the infrared optical conductivity of 
LuFe 2 04 by Fourier- Transform spectroscopy. We per- 
formed reflectivity measurements on a large single crys- 
tal from 7 to 16000 cm -1 (i. c., from about 1 meV to 
2 eV), in the temperature range 10 - 500 K, with the 
radiation field E polarized both ± c and || c. We could 
thus observe in detail the optical phonon spectrum of 
LuFe2C>4, and the striking modifications that it under- 
goes when crossing Tco- In the sub-THz spectral range 
we detected at low T a peak at 10 cm" 1 , which corre- 
sponds approximately to the central mode reported in 
Ref. [l5|. We propose for that spectral feature an alterna- 
tive interpretation which is related to charge order rather 
than to ferroelectricity. 



at frequencies both lower and higher than the measured 
interval. 



RESULTS AND DISCUSSION 
Reflectivity and high-frequency conductivity 

The reflectivity of the LFO crystal is reported in Fig[2] 
for E _L c (a) and E || c (b). In both panels R(u>) shows a 
rich and detailed optical-phonon spectrum, also reported 
on a linear scale in the insets, which changes drastically at 
Tco- In the Elc configuration, at the lowest frequency 
and high T, the non-zero slope of R(<u) reveals a small 
and narrow Drude-like term, thermally activated. Below 
Tco > hi the sub-THz region one instead detects a weak 
but well defined peak at about 10 cm -1 . The mid- and 
the near-infrared ranges are finally occupied by a broad, 
practically T-indcpendent band with a maximum around 
8000 cm" 1 (1 eV). 

This broad peak is better seen in Fig. [3l which 



EXPERIMENT 

The measurements were performed on a large (8x3x2 
mm 3 ) single crystal of LuFe2C>4, grown by optical floating 
zone melting method in a flowing argon atmosphere, [16[ 
and characterized by X-ray diffraction and Laue imag- 
ing at room temperature. The magnetic susceptibility of 
the sample, as measured in zero-fied cooling, is shown 
in Fig. [TJ In addition to the sharp peak at Tjy = 230 
K, also reported previously, [l5| we observe a change of 
slope at about 150 K which may be related to the change 
in the magnetic correlation length observed in neutron 
scattering experiments around 175 K. 

The largest crystal surface contained the axes c and 
a (or b). Its reflectivity R(ui) was measured between 
7 and 16000 cm -1 by Michelson interferometers, with 
the radiation field E polarized either _L c or || c. Be- 
low 300 K the sample was theromoregulated within ± 
2 K in a He-flow cryostat, above 300 K in an evacu- 
ated capsule and in thermal contact with a heater. Both 
sample holders included a mobile hot filament for gold 
evaporation. Indeed, after measuring the intensity I s (oj) 
reflected by the sample, at each temperature gold was 
evaporated in situ, and the intensity Iq{ui) reflected by 
the golden sample was measured for reference. The re- 
flectivity R(uj) = I s ((jj)/Iq(uj) was thus obtained. The 
sub-THz measurements were performed by the same pro- 
cedure but using the coherent syncrotron radiation ex- 
tracted from a bending magnet of the BESSY-II stor- 
age ring when working in its low-a mode. [17J The real 



C* 




part of the optical conductivity was calculated by the 
use of Kramers-Kronig transformations. Standard ex- 
trapolations were applied to R(oj) at all temperatures, 
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FIG. 2: (Color online) Reflectivity of LuFe 2 4 between 10K 
and 500K, as measured with the radiation field polarized 
along the a(&)-axis (above) and along the c-axis (below). The 
phonon region is reported in the inset on a linear frequency 
scale. 
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FIG. 3: (Color online) Optical conductivity of LuFe2C>4 at 10 
K in the whole energy range, as measured with the radiation 
field E _L c or E || c. The fitting curves are reported by dashed 
lines. The zero of the upper curves has been shifted by 200 
r2 _1 cm _1 for sake of clarity 



shows the real part of the optical conductivity, as ex- 
tracted from R(co) by the Kramers-Kronig transforma- 
tions. The Figure reports the data at 10 K, but the 
high-frequency conductivity is virtually independent of 
temperature. The best fit, also plotted in the Figure, re- 
quires a strong oscillator peaked at 7200 cm -1 (0.9 eV) 
for E _L c and at 8400 cm" 1 (1.05 eV) for E || c. The 
former one is in fair agreement with the contribution at 
1.1 eV reported in Rcf. |l4j and assigned to the lowest- 
energy allowed electronic transition, [18[ namely the Fe 2+ 



Fe 3+ charge transfer. This may occur between one of 
the doubly degenerate E levels, which slightly differ in 
energy from each other, and the A-symmctry singlet at 
higher energy. [19| 



The authors of Ref. [1J also report a much weaker con- 
tribution at 0.3 eV, detected in transmission measure- 
ments on thin samples. They assign it to on-site Fe 3+ 
excitations, but it might also be due to the photoioniza- 
tion of the polaronic charges in the CO state. [20| Indeed, 
the mobility of the charges in the CO regime is treated 
in Ref. [lj in terms of small-polaron hopping. In the 
present reflectivity spectra we do observe a weak, lower- 
energy sideband of the main peak, that however the fit 
places (at 10 K) at 4000 cm" 1 (0.5 eV) for E 1 c and at 
4500 cm" 1 (0.55 eV) for E || c. 



Optical conductivity in the phonon region 

The high-temperature, rhombohcdral Rim unit cell 
of LUFC2O4 contains 3 formula units, rotated by 120° 
around c with respect to one another, having the same 



2ll assumed 3 E u and 3 E g (doubly degenerate) opti- 



cal phonons in the ab plane, 3 Ai u and 3 A\ g (non- 
degenerate) optical modes for the c axis, and calculated 
their frequencies (the u modes are infrared active, the q's 
are Raman active). However, a factor group analysis [22( 
leads us to a more complex phonon scenario for LUFC2O4, 
as it is shown in Table I. 



Atoms Wyckoff site 

notation symmetry 



Irreducible 
representation 



Fe, Fe 1(a), 1(b) D id 
Lu,0 2(c) C 3v 

0,0,0 3(d) C 2h 



2E u -\-2 Aiu 
E g + E u +A lg + A 
3E U +A\ U + 2A2v 



TABLE I: Site symmetries of the 7 atoms per formula unit, 
and irreducible representations of the vibrational modes, in 
the rhombohedral unit cell R3m of LuFe2 04. 





Infrared Raman Acoustic 


Silent 




active 


active 






ab plane 


5E U 


E 3 


E u 




c axis 


4A 2u 


Alg 


A 2u 


A lu 



TABLE II: Classification of the phonon modes predicted for 
the ab plane and the c axis, in the R3m high-temperature 
phase of LuFe2 04. 



The resulting modes are classified in Table II with re- 
spect to their polarization and activity. In the high-T 
R3m symmetry, the infrared-active modes are predicted 
to be 5 E u in the ab plane, 4 A2 U along the c axis. At 
low T, the CO reduces the crystal symmetry to the mon- 
oclinic C2/m, [8j which includes one site of symmetry 
C 2/l (l), two Ci(2), two C 2 (2), and two C s (2). Here, the 
number in parenthesis indicates how many atoms each 
site can accommodate. One thus obtains eight possi- 
ble irreducible representations Ti, depending on the site 
occupation. The IYs with i = 1,...,4 lead to 12 one- 
dimensional infrared-active modes, those with i = 5, ..., 8 
to 15 such modes. However, six out of them require that 
a site with inversion symmetry d is occupied. As such 
an atom does not exist even in the high-T phase, they 
can all be excluded. The remaining two irreducible rep- 
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FIG. 4: (Color online) Optical conductivity of LuFe2C>4 be- 
tween 10K and 500K in the phonon region, as measured with 
the radiation field Elc and E II c. 



resentations of the crystal vibrations, in the low-T C2/m 
phase, are 



ri 


= AA U + 8B U + 2A g + 4B g 


r 2 


= AA U + 8B U + 3A g + 3B g 



(1) 

They differ for the Raman spectrum but predict the same 
12 infrared-active, non-degenerate, modes. One thus ob- 
tains a reliable prediction on the effect of the symmetry 
reduction below Tqo on the phonon spectrum. 

This effect is clearly seen in FigHJ which shows the op- 
tical conductivity of LFO measured in the far infrared. 
Below 320 K, for E _L c the broad contribution at 530 
cm -1 splits into two components separated by about 40 
cm -1 ; for E || c all lines shift abruptly to higher fre- 
quencies and new components do appear. Incidentally, 
no phonon line displays Fano-like asymmetries, |23j thus 
confirming the absence of a free-carrier background in 
the phonon energy range. 

In order to obtain a closer comparison with theory, 
the optical conductivity was fit, through the relation 
a = (cli/47t)£2, to the imaginary part e 2 of the Lorentzian 
dielectric function 
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FIG. 5: Optical conductivity at 10K and 500K, for both E 
_L c and E || c, with the fitting curves which provide the 
parameters listed in Table III. The zero of the upper curves 
has been shifted by 200 Sl _1 cm _1 for sake of clarity. 



ei + ie 2 



Aj^TOj 



3 = 1 "TO, 



W 



tJjOJ 



(2) 



Therein, e^ accounts for the high-frequency contribu- 
tions to e(uj) , while Qto > A? ana - 7.7 are t ne central fre- 
quency, strength and damping of the j-th mode, respec- 
tively. The accuracy of the procedure can be evaluated 
in Fig. [5l where the fitting curves are compared with the 
experimental <t(w) at both extreme temperatures. 

The parameters obtained by fitting Eq. [2] to the op- 
tical conductivity of Fig. |4] with the minimum possible 
number of oscillators, both in the rhombohcdral R3m cell 
at 500 K, and in the monoclinic C2/m cell at 10 K, are 
listed in Table III. Therein, a good agreement with the 
factor-group predictions reported in Table II and in Eq. 
[T] is displayed. At high T one obtains a good fit by using 
5 Lorcntzians in the ab plane, 3 along the c axis, to be 
compared with the 5 and 4 infrared modes, respectively, 
predicted in Table II. At low T, in the monoclinic sym- 
metry, the best fit requires 6 modes for E _L c and 6 for 
E || c , in excellent agreement with the 12 infrared active 
phonons predicted in total by Eq. [TJ A shoulder which 
appears at 316 cm -1 both for E _L c and E || c is unlikely 
to be due to a normal mode of the crystal. 

The temperature evolution of the main phonon fre- 
quencies is shown in Fig. [5] for both the ab plane and 
the c axis. Therein one may clearly see the branching at 
Tco of the highest-energy mode of the ab plane and of 
two modes of the c axis. In all these cases, a narrower 
line originates from the high-T absorption at a moder- 
ately higher frequency, while a new mode springs up at 
a much higher frequency 



c 

H 

3 



600- 




M?0 




500- 








400- 


i — • — 


-*-»»l— »- 


— • 


300- 


► — • — 


~* — ••i— — 




200- 


E j_c 

, , . , i , , , , 


a) 




Mode Qto 7 A Q ca i c firo 7 A 



250 

T(K) 



500 



250 500 

T(K) 



FIG. 6: Temperature dependence of the phonon frequencies 
for E _L c (left) and E || c (right) in LuFe 2 4 . 



Optical response in the sub-THz region 

As already mentioned, we have extended our investiga- 
tion well below the phonon region, in the sub-THz range, 
by use of coherent synchrotron radiation. We could thus 
greatly enhance the signal-to-noise ratio on the reflectiv- 
ity, while maintaining a low perturbation on the system. 
The R{u>) thus obtained can be seen in Fig. [2] thanks to 
its logarithmic frequency scale. Therein, at 100 and 10 
K a well defined peak appears at loq = 10 cm -1 . This 
extremely low-frequency feature is then associated with 
the three-dimensional charge order, further reinforced by 
the ferrimagnetic transition at 240 K. The resulting sub- 
THz optical conductivity is reported in Fig. [JJ with the 
feature at ujq again evident at low temperature. In Rcf. 
15l . a peak was detected in the e? extracted from time- 
resolved sub-THz spectroscopy. [15| At 50 K it was re- 
ported at 13 enr 1 (left arrow in Fig. [Jj, in agreement 
within errors with the present ujq. The authors of Rcf. 
15l attributed their peak to a central mode of the fer- 
roelectric LFO phase. They observed a further, broad 
contribution around 40 cm -1 (right arrow in Fig. [Jj, 
which was assigned to a soft mode. In Fig. [Jj at 10 K, 
a second absorption possibly appears between 20 and 30 
cm -1 . 



The present observations may be consistent with those 
of time-resolved spectroscopy, if one considers that in 
polycrystalline samples, as is the case of Ref. [l5j , the 
soft mode stiffens with respect to that observed in single 
crystals. [24j However, if one considers that LuFe2C>4, 
as already mentioned, does not exhibit bulk ferroelec- 
tricity at zero field, there are alternative explanations 
for the present peak at 10 cur 1 (and for the broad 
band at higher frequency, if confirmed by further ex- 
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TABLE III: Parameters obtained by fitting Eq. to the 
phonon spectrum of LuFe2C>4, as measured with the radia- 
tion field polarized in the afe-plane and along the c-axis. The 
calculated values are those reported in Ref. |2U All figures 
are in cm - . 



periments). Indeed, similar features have been observed 
in several systems where the charge order can be de- 
scribed in terms of charge density waves (CDW) and this 
is also the case of LFO. [9( For example, in the manganite 
Lao.25Cao.75Mn03, below the charge ordering tempera- 
ture Tco & narrow peak at 7.5 cm -1 is observed, followed 
by a broad absorption at ~ 30 cm -1 . [25j The former fea- 
ture has been interpreted as a collective excitation (pha- 
son) of the CDW, which displaces from zero to a finite 
frequency when either it is pinned to lattice impurities or 
the charge order is commensurate with the lattice. The 
broader band was ascribed to a combination band of the 
phason with the amplitudon, another excitation which is 
present in a CDW. [25[ According to other authors, the 
highest-frequency feature is instead due to an acoustic 
phonon, which becomes infrared active at zero wavevec- 
tor due to the folding of the Brillouin zone determined 
by the CO. |26[ Further sub-THz experiments are prob- 
ably needed, before reaching a common interpretation of 
the elusive excitations which are detected at those low 
frequencies in the charge-ordered systems. 
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FIG. 7: Optical conductivity for E _!_ c of LuFe204 in the sub- 
THz region at different temperatures. The arrows mark the 
frequencies of the peaks in £2 extracted from a time-domain 
experiment. |15l | 



CONCLUSION 

In conclusion, we have reported here a study of the 
charge-ordered multifcrroic LuFe2C>4 in the infrared and 
the sub-THz range, with special focus on the optical 
phonon region. Therein we have observed the dramatic 
effect of the symmetry reduction - caused by the transi- 
tion at Tco = 320 K between the 2-D and 3-D charge 
order - both on the modes of the ab plane and on those 
of the c axis of the rhombohcdral structure. Below Tco 
some phonon bands split into pairs of narrower lines, 
while new modes do appear. The number of the observed 
phonon lines, both above and below Tco, is in good 
agreement with a factor-group analysis of the respective 
cell symmetries. No evident effects on the phonon spec- 
trum is instead detected when crossing the ferrimagnetic 
transition at 240 K. In the sub-THz region, a very weak 
electric conduction, thermally activated, is detected at 
high temperature, while a peak becomes evident in the 
low-temperature optical conductivity, possibly accompa- 
nied by a broad absorption at higher frequency. For that 
feature, similar to the one reported previously in time- 
domain experiments, we have proposed an interpretation 
related to the charge order rather than to ferroclcctricity, 
and based on a comparison with similar sub-THz spectra 
recently obtained in charge-ordered manganites at low 
temperature. 
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